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Abstract. Irradiating a ferromagnetic material with an ultrashort laser pulse leads 
c/3 , to demagnetization on a femtosecond timescale. We implement Elliott- Yafet type 

spin-flip scattering, mediated by electron-electron and electron-phonon collisions, into 
the framework of a spin-resolved Boltzmann equation. Considering three mutually 
coupled reservoirs, (i) spin-up electrons, (ii) spin-down electrons and (iii) phonons, 
we trace non-equilibrium electron distributions during and after laser excitation. We 
identify the driving force for ultrafast magnetization dynamics as the equilibration 
of temperatures and chemical potentials between the electronic subsystems. This 
*^ ■ principle can be used to easily predict the maximum quenching of magnetization 

upon ultrashort laser irradiation in any material, as we show for the example of 3d- 
| ferromagnetic nickel. 
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Interaction of a femtosecond laser pulse with ferromagnetic metal is known for 
more than a decade to cause ultrafast quenching of the magnetization [TJ. Typical 
demagnetization times extracted from experiments on the 3d-ferromagnetic metals Co, 
Fe and Ni are between 100 fs and 300 fs, depending on the laser fluence and the specific 
sample [21 El El E]- From the theoretical side, intense efforts are being undertaken 
in order to identify the microscopic mechanisms responsible for the observed behavior 
[61 d [U [91 HOl HH [121 UHl US]. Among the considered mechanisms, Elliott- Yafet 
(EY)-type spin-flip scattering [16] is definitely one of the most implemented for the 
modeling of ultrafast demagnetization. The most popular EY-type spin-flip mechanism 
is electron-phonon (el-ph) scattering [21 QUI [131 HH 15] , while recently electron-electron 
(el-el) Coulomb scattering has been considered as well [TT] [TT] . Based on EY-type el- 
ph scattering, Koopmans et al. developed the so-called Microscopic Three- Temperature 
Model (M3TM), which has been shown to possess a great predictive power, explaining 
not only the demagnetization in ferromagnetic transition metals but also the specific 
demagnetization characteristics of gadolinium [2] . Electron-electron Coulomb scattering 
has been implemented in a dynamical model including momentum- and spin- dependent 
carrier scattering [TT], in the following referred to as the el-el Coulomb (EEC) model. 
Since the band structure is taken as input for the calculations, the EEC model can 
be successfully applied to understand the demagnetization dynamics in materials with 
peculiar band structure, as recently shown for the case of half-metallic Heusler alloys 

However, a fundamental question remains still unsolved: What is the ultimate 
"driving force" for ultrafast demagnetization? 




Figure 1. Scheme of the spin-resolved Boltzmann equation: It separates the phononic 
and the electronic system, the latter subdivided into a majority and a minority part. 
Due to the spin-mixing, both electronic systems arc coupled by electron-electron 
and electron-phonon spin-flips. In its essence, the driving force for the ultrafast 
demagnetization is determined by an equilibration of temperatures and chemical 
potentials between the two electronic systems. 
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To answer this question, we introduce a simplified, yet microscopic kinetic model, 
describing non-equilibrium electron dynamics in laser-excited ferromagnetic metals. It 
combines the basic ingredients of the M3TM and the EEC model in the framework 
of a spin-resolved Boltzmann equation. The main idea behind this approach is 
schematically depicted in Fig. [TJ Based on the EEC model, we describe a ferromagnetic 
metal as consisting of three subsystems: the spin-up (majority) electrons, the spin- 
down (minority) electrons, and the phonon system. We consider el-el spin-flips and 
additionally, in the spirit of the M3TM, we include el-ph mediated spin-flips but without 
the need of introducing a separate spin subsystem. In order to enable the implementation 
of complete Boltzmann collision integrals and to ensure a clear interpretation of the 
results it is necessary to make simplifying assumptions, like the free electron gas density 
of states (DOS) for the two spin subsystems. The strength of our approach is that 
it allows to assign a respective temperature and chemical potential to up and down 
electrons, denoted by and //^. We find that the equilibration of these quantities 
is the driving force of ultrafast demagnetization. This equilibration condition can be 
easily applied to real materials with complex DOS to obtain the maximum magnetization 
quenching following laser excitation, as we exemplarily show for nickel. As a further 
result, our approach is capable to clarify the relative role played by el-el and el-ph 
scattering during ultrafast demagnetization. 

In our model, we calculate the time- and energy dependent dynamics of the 
considered subsystems on the basis of a kinetic approach applying complete Boltzmann 
collision integrals [18] . The interactions between electrons (el-el) , electrons and phonons 
(el-ph) as well as the absorption of laser energy (absorp) are described as independent 
collision integrals r co i, respectively. Their sum yields the transient evolutions of the 
electronic distribution f(E) and the phononic distribution g(E). The particular collision 
integrals are generally derived by Fermi's golden rule 

r co i= yl^ol^coil^i)! 2 6{E -E 1 ) , (1) 

all states 

where H co \ indicates the Hamiltonian for the considered collision process, <^ the initial 
state with energy E and ipi the final state with energy E\. Such model was developed 
in Ref. [18] to describe non-equilibrium electron and phonon dynamics during and after 
ultrafast laser excitation for a non-magnetic material. It is capable to account for non- 
equilibrium distribution functions where no real temperatures are defined. 

First, to apply this approach to ferromagnetic metals we include the Stoner model 
into the existing theory. This leads to a separated description of the two spin systems 
(up and down electrons), with a constant energy shift by the exchange energy A ex 
between the corresponding distribution functions. 

Second, we include spin-mixing in our approach. This allows a coupling between 
both electronic reservoirs: Typically, the orthogonality of the up and down states 
guarantees that the matrix element (f | \) vanishes and for a spin diagonal interaction 
operator, there will be no possibility for an electron to change its spin during a collision 
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process. But due to the spin orbit coupling in solids, no pure up or down states are 
defined, but a mixture of both [T3l [T4"] . In this case, an up-electron, for instance, is 
described by the state 

|f) = a |t) + b ||) , (2) 

whereby |6| <C \a\. As the matrix element of the mixed state is larger than zero 
(| (til) I > 0), the spin orbit coupling allows for spin-flip and thus induces an interaction 
between the two spin systems in our model. The collisions including the spin-mixing 
are thus described by 

f co i= -T-\(< 7 o\(<fo\H col \cp 1 )\a 1 )\ 2 5(E -E 1 ) 



h 

all states 

lfa>ki)| a E ^ol^coil^i)! 2 S(E Q -E 1 

all states 

l(o"o|ci)l 2 r co i , 



2 ME„-EA ( 3 ) 



whereby (fi is the spatial and cr, G {t, 1} the (mixed) spin wave function of the electron. 

With all the described ingredients, we can now write down the spin-resolved 
Boltzmann equation 

%r = Yl K/^M)! 2 r ^i 

a ' u ' x (4a) 
+ ^2 \(M\ 2 r ei-ph + ^2 I <^|cr> i 2 r absorp 

cr cr 
§ = El^° r )| 2r Ph-el , (4b) 

with fi, u, cr, A G {t, 1} indicating the spin states. Eq. fll]) consists of three equations 
for f>, and g, respectively. The collision terms r co i can be found in Refs. [HI 119] . 
Each of the distinct equations is coupled to both others, thus the dynamics of all three 
subsystems mutually depend on each other. Since the description allows for spin-flips, 
the number of particles in each electronic reservoir is not conserved. Instead, the sum of 
the transient number of particles of spin-up (nt) and spin-down (?tI) is conserved. The 
transient magnetization M is calculated through 

M(t) = ii B {n\t) - n\t)) = // B An(t) , (5) 

whereby /xg is the Bohr magneton. 

We solve the spin-resolved Boltzmann equation fl3j) for a prototype of a 
ferromagnetic material with the Fermi energy Ep = 8 eV, the effective mass m* = 
1.45 m e and the density of states of a free electron gas. We assume the exchange energy 
A ex = 2 eV. The spin- mixing parameter b 2 ~ 0.03 in Eq. (J2]) is chosen according to 
ab initio calculations for Ni [21 [131 E]- F° r the phononic system we apply the Debye 
model with To = 724 K as Debye temperature. These parameters represent typical 
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Figure 2. Simulation neglecting the contribution of phonons. (a) The internal 
energy U''^ and number of particles of each electronic system, (b) Applying 
Eq. ijSJ) it is possible to define approximated temperatures and chemical potentials 
[V™. The chemical potentials and the corresponding temperatures approach each other 
asymptotically. 



values for real metals. Note that none of them is used as a fit parameter, because it is 
not in the scope of this article to reproduce the experimental data, but to identify the 
driving force of ultrafast demagnetization. In this spirit we chose the density of states of 
a free electron gas because essential relations are given analytically and the solution of 
the spin-resolved Boltzmann equation can thus be easily interpreted. Finally, to clearly 
investigate the processes after optical excitation, we assume a rectangular laser profile 
in time. 

In order to solve the spin-resolved Boltzmann equation fll]) we compute a system 
of approximately 300 strongly coupled integro-differential equations. The solution 
provides the transient non-equilibrium distribution function of each electronic reservoir. 
It resembles the evolution of the electron distribution in Fig. 1 of Ref. [18], but extended 
by the information on the two spin reservoirs. In particular, a net spin-flip is observed 
as a slight particle exchange between both electron reservoirs in the solution of Eq. (J4j) 
during and after irradiation. Figures |2] and [3] were calculated for a 100 fs laser pulse with 
wavelength A = 630 nm and a total absorbed fluence of F — 0.13 mJ/cm 2 . In Fig. |3l 
which will be discussed later in detail, curve (a) shows the transient magnetization 
(normalized to the initial magnetization) obtained with the spin-resolved Boltzmann 
equation (BJ and applying Eq. (jHJ). We obtain a demagnetization time of tm ~ 100 fs 
after laser irradiation. Remagnetization occurs on a picosecond timescale. 

Generally, when a system is excited out of its equilibrium state, it tends back 
to a new equilibrium. An example is the relaxation between electrons and phonons 
after ultrafast laser excitation. Here, the initially heated electrons transfer energy to 
the lattice in a way that both end at a new, higher temperature. The nonequilibrium 
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of temperatures between both systems drives this energy exchange. However, what 
is the driving force of particle exchange between spin-up and spin-down electrons, 
which, according to Eq. (JHJ), is the basis for magnetization dynamics? For such a 
system, allowing for energy and particle exchange, the equilibrium condition refers to 
temperature and chemical potential and reads 

T t = T ; and = ^ + A ex . (6) 

Exploiting the simplifications possible for a free electron gas, the chemical potentials 
can be explicitly expressed by the Sommerfeld expansion: 



fi u (T^) = E F {n u ) 



12 \E F (n^] 1 



(7) 



After optical excitation, the temperatures and will be increased. Let us assume 
that and were already equalized. Then Eq. ([7]) clarifies that in order to equal 
also the chemical potentials ^ and /i^, a change of E F (n^^), mediated by a change of 
rv and n^, is required. Thus, it is the equilibration of chemical potentials of the two 
electronic subsystems, which leads to a particle exchange between both reservoirs and 
thus provides the ultimate driving force for ultrafast magnetization dynamics. 

Our model allows to verify the conclusion of the driving force behind the 
demagnetization process because we can selectively suppress different scattering 
channels in Eq. (j3J). In the following, we discard the phononic influence by setting 
Tci-ph = Tph-ci = 0. The transient spin- resolved density rfl'^ and internal energy 

are depicted in dependence on time in Fig. [2(a) together with the difference An, 
being proportional to the magnetization. Due to the higher electron density of majority 
electrons the amount of energy absorbed by the spin-up electrons is larger than for the 
spin-down electrons. For each moment in time and both spin directions we may define 
a temperature and a chemical potential /r >^ of the current distribution function. 
Both quantities are found through the implicit integral equations of the corresponding 
Fermi distribution f F : 

n u (t) = j fl'\E,^,fi u )BOS(E) dE , (8a) 
U^(t) = j f/{E, T^, //'+) DOS(£) E dE . (8b) 

The resulting temperatures and chemical potentials /i^ are depicted in Fig. [2(b). 
Due to laser excitation, both temperatures increase, while the chemical potentials 
drop. When the laser is turned off (at t = 100 fs) the two spin systems differ in 
temperatures and chemical potentials. Since this difference is small for the assumed case 
of free electron like density of states (second order term in Eq. ([71)), the quenching of 
magnetization is also considerably lower than observed in real materials. Subsequently, 
the reservoirs equilibrate until a same temperature T and chemical potential \i are 




time (ps) 

Figure 3. Magnetization dynamics regarding the different scattering processes 
separately. The calculations were performed by applying the spin-resolved Boltzmann 
equation for a free electron gas. 

reached. According to Eq. (]6]) both constraints have to be satisfied to reach an 
equilibrium between the electronic systems. Fig. |2](b) verifies that the temperatures and 
chemical potentials equilibrate simultaneously. Comparing with Fig. E](a), we find that 
also the transient magnetization, being proportional to An, reaches its asymptotical 
value on the same timescale. In its essence, the particle exchange affected by the 
difference in chemical potentials provides the driving force for ultrafast magnetization. 

In order to estimate the effect for more realistic systems, we choose nickel and solve 
the equilibrium condition fl6]) applying the implicit equations (JEJ) with the real density 
of states for Ni [20]. Note that this can be done without explicitly solving the spin- 
resolved Boltzmann equation. In case of starting at room temperature and ending up 
at an electron temperature of T e = = = 3000 K a maximum quenching of 30% 
is found. Such transient electron temperatures are typical for ultrafast laser excitation 
of metals below melting threshold [181 El]- Thus, applying the main result of our 
simplified model to realistic systems leads to quantitative agreement with experiments, 
demonstrating the solidity of our findings. 

To show the great potential of the spin-resolved Boltzmann equation, we now 
exemplarily compare the relative effects of el-el spin-flips and el-ph spin-flips on the 
demagnetization process. We again restrict ourselves to the case of a free electron 
gas. Figure |3] represents the calculated magnetization dynamics obtained with the 
spin-resolved Boltzmann equation ((4]), when different collision processes are regarded 
separately. The red solid line (d) corresponds to the red solid line in Fig. [2](a), 
determining the magnetization dynamics when the influence of the phonon bath is 
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completely discarded. Including the cooling by phonons but disregarding el-ph spin- 
flips (green dotted line (c)), we see a slight recovery of magnetization. 

Thus, phonon cooling is essential for the remagnetization process. However, 
phonons also strengthen the magnetization dynamics: If we additionally allow 
for phonon mediated spin-flips we find the black dash-dotted curve (a). The 
demagnetization occurs faster than in the case where only electron collisions mediate 
spin-flips, while the remagnetization is strongly accelerated by phonon mediated spin- 
flips. Hence, phonons can act like a catalyst for the ultrafast magnetization process. 
The blue dashed line (b) is calculated excluding electron-electron-collision mediated 
spin-flips. In this case, the quenching of magnetization will be less. Comparing the 
simulations shown in Fig. [3] we conclude that for our system electron mediated spin- 
flips are important mainly during the demagnetization phase, while phonon mediated 
spin-flips dominate the remagnetization phase and increase the effectiveness of the 
demagnet izat ion . 

In conclusion, we described the ultrafast demagnetization with the spin-resolved 
Boltzmann equation. We identified the driving force of the demagnetization process 
in the equilibration of the temperatures and chemical potentials of the up and down 
electrons. As shown for the example of Ni this equilibrium condition provides a 
possibility to easily estimate the maximum quenching for any ferromagnetic material. 
Our approach is also capable to investigate the role of different scattering processes 
during ultrafast demagnetization dynamics. 
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